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The matching problem between functional shapes via a T?V-penalty 

term: a T-convergence result 

B. Charlier * G. Nardi, A. Trouve ** 


Abstract 

In this paper we study a variant of the matching model between functional shapes introduced 
in [B]. Such a model allows to compare surfaces equipped with a signal and the matching energy 
is defined by the L 2 -norm of the signal on the surface and a varifold-type attachment term. 

In this work we study the problem with fixed geometry which means that we optimize the 
initial signal (supported on the initial surface) with respect to a target signal supported on a 
different surface. In particular, we consider a BV or U 1 -penalty for the signal instead of its 
L 2 -norm. Several numerical examples are shown in order to prove that the BV-penalty improves 
the quality of the matching. Moreover, we prove a T-convergence result for the discrete matching 
energy towards the continuous-one. 


1 Introduction 

Previous works. In jl>], the authors introduce a new framework to study the matching problem 
for surfaces with boundary equipped with a signal (functional shapes or fshapes) called fshapes. New 
developments in non-invasive acquisition techniques sue as Magnetic Resonance Imaging (MRI) or 
Optical Coherence Tomography (OCT) allow to get geometric-functional data for several diseases 
(e.g., cortical thickness in the study of Alzheimer’s disease or thickness of retina layers for the evolu¬ 
tion of glaucoma). Then, new methods for the registration of fshapes are needed in medical imaging 
in order to describe anatomical variability and produce statistical estimate for several diseases. 

The energy matching considered in |6j takes into account the L 2 -norm of the signal on the surface 
and a varifold-type attachment term. Tools from geometric measure theory (currents, varifolds) are 
increasingly used in image processing m They actually allow to define a distance between two 
geometric objects in the setting of the measure theory. Moreover, because of their properties of 
compactness, they represent the good setting in order to define energies to minimize. However, such 
distances are defined in a weak sense and their computation can be very difficult. For this reason the 
model in |Bj is developed in the setting of the Reproducing Kernel Hilbert Spaces (RKHS), which 
allows to easily compute a varifold distance. We refer to [IS] and to Section [2] for more details. 

The main difficulty of such a model arises from the fact that the position and the signal are 
minimized simultaneously. Now, the modification of one of these variables (position or signal) can be 
sufficient to get a local optimal configuration. Numerical simulations show that the gradient descent 
can be trapped in local minima where the signal has been strongly modified and the geometry has 
not. 
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This phenomenon is essentially due to two facts: the choice of the penalty term for the signal and 
the choice of Gaussian kernels to estimate the varifold-norm (see Example |3.4[ ). The main object of 
this paper is to show that a -BF-penalty term can strongly improve the matching results. 

Contributions. In this work we consider the “functional” matching problem between (X, /) and 
a target fshape ( Y,g ). We consider the following three energies: 

E(f) = R(f,X) + V((X,f),(Y,g))) 

where R(f, X) can be defined as 

R(f,X) = l\\f\\h(X) (L 2 - model), 

R(f,X) = l\\f\\ 2 Hl{x) (H 1 - model), 

R(f,X) = \\f\\ BV{x) (BV-model). 

and the attachment term V is defined by using the varifold theory (see Section [ 2 ]). The L 2 -model 
represents the model introduced in (6|. 

We are interested here in the minimization problem with fixed geometry. This means that the 
optimization is made only with respect to the signal. In other words, the optimal fshape is supported 
on the initial surface X. We prove in particular an existence result for the optimal signal in the case 
of the BV and I5I 1 -model. The existence result for the L 2 -model is already proved in |BJ theorem 6. 

Afterwards, we define a discrete version of the problem by approximating the surface X by a 
sequence of triangulations. The continuous problem can be approximated by a sequence of discrete 
problems defined on some triangulations whose triangles’s diameter goes to zero. Roughly speaking, 
smaller is the diameter of the triangles higher is the number of the vertices and, when the diam¬ 
eter goes to zero, the triangulation converges to the initial surface (with respect to the Hausdorff 
distance). 

Concerning the definition of the discrete problems, the main issue is represented by the choice of 
the admissible triangulations. In fact, as X has a boundary, the triangulation and the surface need 
not be one-to-one. In this paper, we decided to work with the class of triangulations that cover the 
surface and whose surpassing part has small area. 

Moreover, we must guarantee that the discrete solution is a ’’good” approximation of the contin¬ 
uous local minimum. This can be proved by the T-convergence theory, that is a natural notion of 
convergence of functionals allowing to justify the passage from discrete to continuous problems. In 
particular, in the case of minimization problems, the T-convergence guarantees also the convergence 
(in some sense) of the discrete minima towards the continuous-one. 

We prove a T-convergence result showing that the minimum of the discrete problem is close to 
the minimum of the continuous problem if the diameter of the triangles is small. The main issue 
to get such a result arises from the fact that the discrete problems and the continuous problem are 
not defined on the same geometric support. Now, in order to get the T-convergence result we need 
some hypothesis on the triangulations which guarantee the convergence of the areas. There is in 
fact a famous example (see Schwartz’s lantern example in [23! p. 354) proving that the area of the 
triangulations described above needs not to converge to the area of X. 

Following US!, we point out the suitable properties a sequence of triangulations should verify to 
guarantee the convergence of the areas. To this end the angle between the normal to the triangu¬ 
lations and the respective (in the sense of the projection on the surface) normal to the surface has 
to go to zero. For every penalty (L 2 , H 1 , BV), previous condition on the triangulations allows to 
prove the T-convergence result for the respective energy. 
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This kind of condition is involved in the numerical study of several problems defined on surfaces 
(e.g., Laplace-Beltrami operator). Then our result can be useful behind the matching model for 
fshapes. 

Finally, we show some numerical examples to compare the different models. These examples 
actually point out that the RF-model strongly improves the matching result. 

Relationship with [6]. The problem studied in this work corresponds to the model presented in 
|6J with different signal-penalty term and fixed geometry. In this paper, we point out the influence 
of the penalty term on the quality of the matching and the properties needed to get a F-convergence 
result. The present work represents a important step in order to improve the general model presented 
in p. 

Moreover, studying the problem with fixed geometry does not make the problem trivial. In 
fact, we perform the matching between two different fshapes that are not supported on the same 
surface. Such a geometric difference is taken into account by the varifold attachment term. Then, 
the influence of the penalty term and its interaction with the varifold term is a interesting problem 
to study. 

In particular, we address the problem of the existence of a solution, the discretization of the 
problem and a related T-convergence result. These problems represent are needed to get a similar 
result in the general case of fshape (i.e., the surface A can vary). The effect of the penalty term on 
the geometric optimization will be addressed in a future work. 

Structure of the paper. In Section [2] we remind the definition of varifold and the RKHS frame¬ 
work used in [61 in order to estimate the varifold norm. In Section [3] we introduce the three models 
(L 2 , H 1 , BV ) studied in this paper and prove the existence result for the BV and iT-models. In 
Section [4] we define the class of admissible triangulations and a condition on the angles between the 
respective normals in order to get the convergence of the area. In Section [5] we define, in the setting 
of the finite elements method, a discrete version of the energies introduced in Section [3} Section [6] is 
devoted to the definition of the discrete problems and the results of T-convergence. Section [7] shows 
several numerical examples in order to compare the different models. Finally, in the Appendix, 
we recall the definition of the BV functions on manifolds and prove that such functions may be 
approximated by C 1 functions. 

2 The functional varifold setting 

In the following we assume that : 

Hypothesis 1. A" is a 2-surface of M 3 with boundary. X is smooth (C 2 at least) oriented, connected 
and compact. 

The boundary of X is denoted by dX and the interior is denoted Ay. The smoothness condition 
in Hypothesis [l] means that Ao is a two dimension smooth manifold and dX is a one dimensional 
manifold of the same regularity. In this section we remind the main tools needed to set the model 
introduced in |6j. 

Definition 2.1 (fshape). We define a fshape as a couple (AT,/) where X is a surface verifying 
Hypothesis [l] and / : X —)■ M is a signal defined on X. 

We denote by L p ( A), W l,p { A), and BV(X) the usual Lebesgue, Sobolev and bounded varia¬ 
tion functional spaces defined on the surface A (see Appendix [9] for details). In order to define a 
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suitable attachment term we remind the definition of functional varifold. The space of continuously 
differential functions defined on a norm space E and vanishing at infinity is denoted Cq(E). 

Definition 2.2 (fvarifolds) . A 2 -dimensional functional varifold (fvarifold) is any linear form be¬ 
longing to (Co (R 3 x G(3, 2) x R)) 7 , the dual of the space of real valued continuous functions vanishing 
at infinity. We denote by G(3, 2) the Grassmannian of the non-oriented 2-dinrensional linear sub¬ 
spaces of R 3 . This space (diffeomorphic to the projective space P 1 (M 3 )) has a structure of smooth 
manifold and can be equipped with a smooth distance 

d(S,T) = (2(1-(n s ,n T ) 2 )) 1/2 

where ns and n t are respectively unit normal vectors to S and T in G(3, 2). 

In particular, we can associate a fvarifold with a fshape (X, f) by setting 

/*(*,/)(¥>) = / V(x, T x X, f(x))dU 2 (x) Vip 6 Cq(M 3 x G(3, 2) x R), 

J X 

whrere (I'M. 2 is the 2-dinrensional Hausdorff (or volume) measure. Remark that, according to Theorem 
1 p. 49 in m, such a functional is represented by a measure on Cq(M 3 x G( 3,2) xR). Next lemma 
is very usefull in the following. It can be easily shown by the previous definition of fvarifold. 

Lemma 2.3. Let { Xh} be a sequence of surfaces such that 7i 2 (Xh) —> 0. Then l-i(x h ,f h ) converges 
to the null varifold for every sequence {//,} of signals. 

We refer to [22] for more properties of varifolds and to [6] for more properties of functional 
varifolds. 

In |6), the authors evaluate the varifold norm by an embedding into a RKHS generated by 
Gaussian kernels. We refer to [24] for an introduction to RKHS theory and to Section 3.2.2 in [6] 
for a general discussion on the construction of a suitable RKHS for our framework. 

Following [5j, we consider W as the RKHS associated with the positive kernel 

k((x,T ly f 1 ),(y,T 2 , f 2 )) = k e (xi,x 2 ) ®k t (T l ,T 2 ) ®k f (f 1 ,f 2 ) 

where k e ,kt,k v are some positive kernels on R 3 , the Grassmannian, and M, respectively. In our 
applications we use the following smooth Gaussian kernels 

|| Z1 -x 2 || 2 2(l-<» Tl . rar ,,) 2 ) I/1-/2I 2 

k e (x,y) = e ^ , k t (Ti,T 2 ) = e , kf(a,b) = e (2.1) 

where a e ,at,<Jf are three positive constants and nx k represents the unit normal vector to T*,. Remark 
that, we define a kernel on the Grassmannian by identifying the linear space (plane) with its normal 
space (line). 

Then, W is continuously embedded in Cq(M 3 x G(3,2) x R) (see Proposition 3 in 0), and we 
can define a continuous map 


i* : (C^R 3 x <7(3,2) x R))' ->■ W'. 
In particular, for every n G (Cq(R 3 x G( 3, 2) x R))' we have 


b\\w' ^sup {n(if) : ||^||c 0 i(K3xG ( 3,2)xK) < !} 


(2.2) 
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where, by an abuse of notation, we denote i*(y) by g. Moreover, as W and W' are Hilbert spaces, 
every functional g G W' can be represented by an element of W via the isometry K\y : W' —> W 
given by 

(K w g, ip) w = g(tp) Mip G W , G W' 
so that, for every fvarifold, we get 

Mw> = \\KwlA\w V/i G (^(M 3 x G{ 3, 2 ) x R))' . 

In particular, the kernel represents every Dirac delta functional, which means that 

K w$( x i,Vi,/i) : {x,V,f)h^k((x 1 ,V L ,fi),(x,VJ)) 
for every (xi, Vi, /i) G M 3 x G(3, 2) x M. Then 


( s (x 1 ,v 1 ,f 1 )^(x 2 ,V2,f2))w' = k e (xi,x 2 )kt(Vi, Vfikftfi, /2), (2.3) 

for every (xi, Vi, f\), (x2, V2, /2) G M 3 x G(3, 2) x M, so that, for every couple of fshape ( X , /), (Y, < 7 ) 
we get 


(M( xj),k-(Y, g ))w' = [_ f k e {x,y)kt{T x X,TyY)kf{f(x),g(y))d'H 2 (x)dn 2 {y). (2.4) 
Moreover, ||/i(.Y,/) — l- l (Y,g) IIW 7 ' is a distance between fshapes which can be easily computed using 
order to define an attachment term between fshapes. 

3 The matching energy: an existence result 

In this section we define the matching energy between two fshapes and we prove an existence 
result for the optimal solution. 

Let be (Y, g ) a target fshape. We consider the following three energies defined for a generic fshape 

E(f) = R(f,X ) + ^I \VKx,f) - tHY,g)\ \w' (3-1) 

where R(f, X) can be defined as 

R(f,X) = l\\f\\h(x) (L 2 - model), 

R(f,X) = k\\f\\ 2 Hl{x) (H 1 - model), 

R(f,X) = ll/ll bv(x) (BV-model) . 

In particular, we aim to solve the minimization problem 

mf E(f) 

where the infimum is taken on the respective space (L 2 (Y),H 1 (Y) or BV(X)). In this section we 
state and prove an existence result for such a problem. In the case of the HY-model and iT-model 
the proofs are very similar so we prove the result for the HY-model and we indicate how it can be 
adapted to the iT-case. 


(2.4). Such a framework (fvarifolds and RKHS dual norm) is used in [ 6 ] and in the present work in 
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We remind that we optimize only with respect to the signal which implies that the initial config¬ 
uration and the optimal one have the same geometric support. However, the geometry is taken into 
account in the attachment term. 

As pointed out in the introduction such a problem is not trivial and strongly depends on the 
choice of the penalty term R. In fact, as we work with Gaussian kernels, an oscillating signal can 
be more economic in terms of energy E (see Example 3.4). This is due essentially to the fact that 
e -\a-b\ j g sma p if |q — 5| is large. The choice of a penalty term which takes into account the gradient 
of the signal prevents this kind of optimal configurations. 

A sequence {nh)h of W is said to weak-* converge to // £ W if 


lim i u h (ip) = n(<p), 

h —>0 


VV e C^(R 3 x G(3,2) x 


In that case we write ^ g,. The following Lemma links the pointwise convergence of a signal 

and the weak-* convergence of fvarifolds. It can be easily proved using the definition of varifold and 


applying the dominated convergence to (2.4). 

Lemma 3.1. If fh -7 / a.e. on X then /i(xj h ) —*■ x,f )• Moreover, 

\\V(x,f)\\w' = ^™\\^x,f h )\\w' ■ 


3.1 The L 2 -model 

The L 2 -model is studied in [B] part I. The authors consider the matching energy 

E (f) = 'y ll/lli 2 (x,K) + Y~\\^{x,f) - V(Y,g)\\w' 
where 7/,7w are two positive constants. They prove the following general result : 


(3.2) 


Theorem 3.2 (Proposition 7 in [6]). Let X,Y be two finite volume bounded 2-rectifiable subsets 
of M 3 . Let us assume that W is continuously embedded in (7o(M 3 x G(3,2) x M) and g G L 2 (Y). 

If the ratio 7 //tw is large enough, then there exists at least one solution to the minimization 
problem 

inf E(f) 
feL 2 (x) 

and every minimizer belongs to L°°(X). Moreover, if X is a C p -surface and W ^7 C'q*(M 3 x G(3, 2) x 
M) with m. ^ max{p, 2}, then every minimizer belongs to C P ~ 1 (X). 

Finally, there exists a constant C > 0 (independent of X and Y) such that, every minimizer 
verifies 

\\fh°°(x) < C^(H 2 (X) + U 2 (Y )). (3.3) 

7 / 

We point out that, for every minimizing sequence (fh)h, a bound on E(fh) guarantees only the 
L 2 -weak compactness for the signals, which is not enough to get the semi-continuity of the varifold 
term (see Lemma 3.1). To avoid this problem the authors relax the problem to the class defined 
as it follows: 

Definition 3.3. A4 x is the class of the Borel finite measures v on M 3 x G{ 3, 2)xR such that 
J tp(x,V)du(x,V,f) = f (p( x ^ T xX)dn 2 ( x ) \/ip g C c (M 3 x G(3,2)). 

Note that £ Xi x for every fshape (A,/). 
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In this framework, they relax the energy E to the following energy 

E:M X ^R, %) = | J\f\ 2 du+^\\v-^ g) \\ 2 wl . (3.4) 

It holds that 

£(/) = %!,/))• (3-5) 

This allows one to get compactness in the setting of measures by using the Prokhorov’s theorem. 
Moreover, it can be shown that the minimizing measure v* of E is actually associated with a fshape, 
so that u* = fitxj*) f° r some f* £ L 2 (X). The proof relies on the Implicit Function Theorem which 
needs the hypothesis on 7 // 7 W. We refer to Proposition 7 and Lemma 2 in [5] for more details. 

Such a model has two main issues. Firstly the existence result depends on the weights used to 
define the energy. Moreover, the L 2 penalty does not prevent from some oscillating configurations 
for the optimal signal. Next example, which is similar to simulations in [Bj Section 9.1, points out 
these facts. 

Example 3.4. We consider the two fshapes shown in Figure [l] two bended smooth stars both lying 
on the cylinder {(x,y,z) £ M 3 |x 2 + y 2 = 1}. The source fshape is a fat star with constant signal 
equals to zero composed of 22880 vertices (45171 triangles). The target fshape is a thin star with 
a signal taken from “Lena” picture widely use in the image processing community (the signal has 
been normalize to range between 0 and 1) composed of 4056 vertices (7744 triangles). 





(a) Source (. X , 0) 


(b) Target (Y,g) 


(c) Relative position (mesh view) 


Figure 1: The initial fshape is a fat smooth star la and the target fshape (Y,g) is a thin star with 


Lena picture as grey level signal [Tb] (black is 0 and white is 1). Remark the pixel effect due to the 
P0 discretization (see Section [5]). The blue cylinder is just drawn to help to see the perspective. 


The two fshapes are “overlapping” although we emphasis that they do not contain the same 
number vertices and do not share any vertices (see Figure lc). In order to minimize the varifold 
norm appearing in the matching energy (3.2), the fshape source has to be carried with a Lena signal 
on the overlapping part of the two fshapes. On the remaining part of the source, the optimal signal 
is determined by the penalty term and it should ideally remain to zero (no information added). 
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We use the fshapesTk software [7] to minimize a discrete version of (3.2). Thorough details 
about the discrete energy is given below. We choose 7 w = 1 since the balance between terms is only 
determine by the ratio jf/jw- We recognize the same behaviour as in the example of | 6 j Section 
9.1: the optimal signal is strongly depending on 7 f and is oscillating on the non overlapping part 
(the Moire pattern visible at the end of the branches on Figures 2b 2c and 2d). The intensity of 
the oscillations as well as the contrast of the optimal signal depend on the value of the parameter 


7 j as shown by Figures 2b and 2c As discussed in [ 6 j, a possible trick improving the quality of the 
optimal signal is to smooth the gradient used along the gradient descent (see Figure 2d). 

As pointed out at the beginning of this section, this essentially depends on the behavior of 
Gaussian kernels. As shown further, this kind of configurations is avoided by the RF-model. 


3.2 The 517-model 

We consider the following energy 


1 , 


E(f) — ll/llW(X) + xll H(X,f) - P{Y,g)\\w' 


(3.6) 


where (Y, g) is a target fshape. The H^-model can be studied in a similar way (see Remark 3.6). 
We are interested in solving the following problem 


inf E(f) . 
feBV(x) 


(3.7) 


We prove now the existence of a solution to problem (3.7). 


Theorem 3.5. Problem (3.7) admits at least one solution. 


Proof. Let (fh)h be a sequence of RF( X) minimizing for E. We can suppose that ||/ft,||sy(x) is 
uniformly bounded and, by the compactness theorem for RF-functions, {fh)h converges (up to a 
subsequence) to some / E BV(X) with respect to weak-* topology. This implies in particular that 
(up to a subsequence) fh —> f a.e. on X. The result ensues by remarking that the fvarifold norm is 


lower semicontinuous with respect to the a.e.-convergence (Lemma 3.1), and that the RF-norm is 
lower semicontinuous with respect to the weak-* topology. □ 

Remark 3.6 (The R 1 -model). We point out that we get the same result if we consider signals 
belonging to R 1 (A) instead of BV(X). This follows from the fact that the unity ball of R 1 (A) is 
compact with respect to the weak topology. Moreover, H 1 (X) is compactly embedded in R 2 (A), 
which implies in particular (up to a subsequence) the a.e.-convergence of the minimizing sequence. 


4 Surfaces and triangulations 

In this section we remind the main properties of triangulations and describe how to compare a 
signal defined on a surface with a signal on a related triangulation. 

4.1 Triangulations of a surface 

We start by stating some general notions about triangulations related to the surface X. We 
remind that X verifies Hypothesis [lj Following [HEDED] we give the following definition of trian¬ 
gulation: 

Definition 4.1 (Triangulations). A triangulation T is a two-dimensional manifold (with bound¬ 
ary) consisting of a finite set A 7- of topological triangles such that: 
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(a) The target (Y, g) 



(b) Optimal fshape ( X , /*) with 7 = 0.1. 



(c) Optimal fshape ( X,f *) with 7 = 1 (d) Optimal fshape (X,/*) with 7 = 0.01 and with 

a smoothed gradient descent. 


Figure 2: Results for the L 2 penalty: en face view of the target (Figure 2a) and of the optimal 
signals f* defined on X (Figures 2b 2c and & The optimal signal f* is obtained by minimization 
of the discrete version of the functional (3.2) with various parameters. 
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1. any point jjgT lies in at least one triangle T E A 7 -; 

2. each point has a neighbourhood that intersects only finitely many triangles T E A 7 

3. the intersection of any two non-identical triangles T, T' E A 7 - is either empty, or consists of a 
common vertex, or of a simple arc that is an edge of each of the two triangles. 

In the following we denote by T a generic triangle of A 7 - and by dT the boundary of the manifold 
T. We call vertices of T the set of vertices of the triangles T E A 7 -. Finally, we define 

diarn^ = max {diam(T)}, (4.1) 

i.e. diarn^ is the supremum of all the diameters of the triangles composing T. 

Definition 4.2 (C-Regular triangulation). We denote by TV/ the set of triangulations T such 
that there exists a constant C such that 

— VTe A t (4.2) 

Pt 

where h t is the diameter of T and px is the diameter of the sphere inscribed in T. 

We define the distance function as 

Vr E K 3 , dx(x) = d(x,X) = inf \x — y\ . 

y&X 

For every x E M 3 , we call (if it exists) projection of x on X every point tvx(x) E X such that 
dx(x) = \x- n x (x)\. 

Definition 4.3 (Tubular neighborhood). Let X be a surface satisfying Hypothesis [lj We denote 
by U r (X) a subset of M 3 of the form 

U r = { x E M 3 | dx{x) < r } 

such that every point x E U r admits a unique projection ttx{x) E X. We refer to m for the proof 
of the existence of such a tubular neighborhood for some r > 0. 




Figure 3: A surface X sketches by a blue curve and its tubular neighborhood. 
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We remark that the projection is not injective on the set U r (X ) but only on the following subset 

N r (X ) = {x + tnx(x) : t €.] — r, r[, x € X} (4.3) 

where we denote by nx(x) the unit normal vector to X at x (see Figure pi). Remark that normal 
vector can be defined on dX by continuity. We also point out that dx € C 2 (J\f r (X)). If x € M r (X) 
then its projection on the surface is defined by the following equation: 


x = 7 r x (x) + d x (x)n x (7rx(x )). 


(4.4) 


As pointed out in the introduction, because of the curvature of dX, we can not assume that 
a triangulation, obtained by the sampling of some points of A, is in bijection with the surface. 
Then, in order to correctly define a discrete setting, we introduce a particular class of triangulations. 
Before stating the definition we need the following notations : let r > 0 and for every triangulation 
T belonging to the tubular neighborhood U r (X) of a surface X we denote 

T n = T n Nr{X) and T out = T n N r (X ) c . 


Definition 4.4 (ft,- admissible triangulations for a surface). Let h > 0. We say that a trian¬ 
gulation T is ft,-admissible for the surface X if the following properties hold 

(*) T lies in Uh{X ); 


(ii) T lies in Nh(X +ri ) for some rj > 0, where X +v is an extension of X defined in (9.14); 
(in) T m and X are in one-to-one correspondence through nx', 

(■ iv ) / H 2 (T ont ) = 0(h). 


(v ) diam- 7 - = 0(h ) defined in (4.2) 



(a) ft.i-admissible triangulation 7^, for X 



(b) ft- 2 -admissible triangulation Th 2 for X 


Figure 4: A surface X with boundary (solid grey) and two admissible triangulations (black lines) 
with h± > li 2 - The area of 7/° 2 ut (i.e. the upper part of triangulation exceeding X) is smaller than 

■T-OUt 

'h 2 ■ 


Previous definition guarantees that every /i-admissible triangulation T can be decomposed in 
two parts: T = T m U T out . The first-one, T m , lies in J\fh(X) and is in bijection with X. We can 
define a partition of X by considering the family {T?;} defined by 

T i = ir x (T i nAr h (X)) MTi G T. 
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The second-one, T out , lies outside A4(A) and condition (iv) guarantees that its area is small (see 
Figure [4]). 

Remark 4.5. In general, an inscribed triangulation (in the sens of m) is not in bijection with a 
surface with smooth boundary through the normal projection. The bijectivity can fail close to dX 
because of the curvature of dX as depicted in Figure 5a Locally, the normal projection of dX on 


a hyper-plane needs not to be a line. In that case, it is not possible to project dX on an edge of a 
triangle. 



(a) Inscribed triangulation 7i for X (b) An ii-admissible triangulation 7^ for X 


Figure 5: The surface X is a bended smooth star (solid grey) as in Figure lb Two triangulations 


(black lines) are pictured. Figure [5a] the triangulation 71 is not in one-to-one correspondence with X 
through the projection map (for instance the part of the smooth star in red exceed the triangulation). 
Figure 5b the subset T™ of 71 is in one-to-one correspondance with X. 


Remark 4.6 (Real data). We point out that the previous definition is slightly different from 
the usual notion of triangulation associated with a manifold. We actually do not suppose that the 
vertices of the T belong to the manifold X. See Figure 5b 


Such a generalization corresponds better to the features of the experimental images. In fact, the 
original images (OCT, functional MRI) are modified (segmentation, deblurring, denoising) in order 
to improve their quality and to be used for the numerical tests. Then, the real data that we use 
represent an approximation of the real objects we want to compare, and, moreover, can contain some 
noise. By considering a set of triangulations whose vertices do not belong to X, the results set in 
this paper are more adequate to such a situation. 


4.2 Convergence of areas 

We remind that the Hausdorff distance between two surfaces X , Y C M 3 is defined as 


dy.(X, Y ) = max ] sup dy{x), sup d\(y) 

1 x&X y&Y J 


If T is an /i-admissible triangulation of X for some h < 0 then property (z) of Definition 4.4 implies 
that du(X,T) < h. Unfortunately, the convergence with respect to Hausdorff distance does not 
imply the convergence of the areas. A counterexample is given by the Schwartz’s lantern (see {23] 
page 354). In the following we remind a result from |I2] giving a sufficient condition in order to have 
the convergence of the areas. 


Definition 4.7 (Angle between normals). Let h > 0. Assume that T is an Ii-admissible trian¬ 
gulation for X. For every x G T m we define the angle a x as follows 
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• if x belongs to the interior of some triangle, then a x is the angle belonging to [0,7r/2] between 
the two normals nx{j rx0*0) an d n T( x ); 

• if x belongs to an edge of a triangle, then a x is the biggest angle (belonging to [0,7r/2]) between 
n x(itx(x)) and the normals of the triangles containing x. 


In the following we set 


^max — Slip OL x • 
xeT in 


Lemma 4.8. Let h > 0 and T be an /r-admissible triangulation for X. We have 


\n 2 (X) - H 2 (T ia )\ = 0(a L x + d n (x, T )) 


Proof. For every x £ X we consider on the tangent space T x X the basis B(x) = (e 1 (x), e 2 (x)} given 
by the two principal directions. We denote by ki(x) and k 2 (x) the principal curvatures of X at x. 
Similarly, for every x £ X we can consider the basis B(x) = (e 1 (x), e 2 (x), nx(i)} for M 3 . 

Consider the differential Dttx ■ M 3 —>• T nx ^X of 7 tx- Note that, for every x £ A fh(X), we have 
Dnx(x)(v) = 0 for every variation v in the direction nx(^x(x)) orthogonal to T nx / x \X. So, we 
should just consider the tangential variations in order to calculate the Jacobian of the projection. 
In [T9] it is proved that for any x £ U^f X) and v parallel to T nx ^X we have 

Dir x (x)(v) = [I\ T ~ - £ x d x (x)Dn x {nx(x)) ) v 

V ) 


as q Ly i dy_(X, T) —> 0 when h —> 0, where a max is introduced in Definition 4.7 


And the matrix of Dirx(x) written with respect to the basis B(nx(x)) and B(ttx(x)) is 




Dtt X {x) = ( 1 + d x{x)£xKl{.^x{.x)) 


0 

_ 1 _ 

l+d X (x)£ x K2(ltx{x)) 


0 ^ 

0 , 


where e x = ( n x (ir x (x))) £ {-1,+1}. 

Let now A C T m be a subset of a triangle of T- This implies that the Jacobian of the projection 
on X restricted to A is given by 


det(L> rA 7Tx)(x) 


cos a x 

(1 + d x {x)£ x Ki('K X {x)))(l + d x (x)£ x K 2 (7rx(x))) 


V x £ A. 


We have cosa x = 1 + 0(a. f. AV ). Moreover, as the principal curvatures are uniformly bounded 
and dx(x) = 0{du{X,T)), we get 


det(D u ir x )(x) = (1 + 0(a 2 max ))(l + 0(d n (X,T))) = 1 + 0« ax + MX, T)) (4.5) 

as a ma xjdn{X,T) 0. The proposition follows by changing the variable in the formula of the 
area. □ 


In view of Lemma |4.8| let us introduce the following set of Assumptions: 

Hypothesis 2. Let {Th)h be a sequence of triangulations indexed by a parameter h > 0 decreasing 
to 0 such that 


(*) for any h > 0, the triangulation Th is /r-admissible for X; 
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(ii) the sequence a 


h 

max 


= 0(h) where ax is the angle of Definition 


4.7 


for the triangulation If. 


Remark 4.9. The parameter h > 0 is both the index of the triangulation and a measure of the 
precision of the triangulation. The idea is to consider decreasing h yielding to a more and more 
accurate. 


This implies in particular the following result: 

Proposition 4.10 (Convergence of the area). Let X be a surface satisfying Hypothesis [I] and 

(Th)h a sequence of triangulation satisfying Hypothesised Then we have lim PL 2 {Th ) = PL 2 {X). 

/i—>o 


Proof. The proof follows from Lemma 4.8 and (iv) of Definition 4.4 


□ 


Proposition 4.10 is a generalization of Corollary 5 of m to surfaces with smooth boundaries. 


4.3 From the triangulation to the surface 

In this section we define how to carry a signal from the triangulation to the surface. 

Definition 4.11 (Projection). For every function / defined on an admissible triangulation for X 
we define the projection of / onto X by 

f e :X^R, f e (7T X (x)) = f{x) x G V n . 

We point out that the function y carries on X the signal defined on T m . 

Proposition 4.12. Let h > 0 and T be and an admissible /i-triangulation for X. Then, for every 
/ G W 1,0 °(T, M), we have 


\\f e \\LP(X) ~ ll/lliP(T in ) ^( a max "h 7 "")), 
l|vVii W ) = l|v ri 7llLr (Tte ) + o(aLx + dw(x,r)). 

for every p G [1, oo], as d ni x , T) -4 0. 

Proof. The first equality is proved by performing the change of variables y = 7r x(x) and using 
Lemma |4.8[ The second relationship is proved by applying the chain rule. A rigorous proof is given 
in tm Lemma 3. □ 


5 Discretization 


5.1 Notations 


Let T be a triangulation in the sens of Definition 4.1 We denote by N v , N e , and Nt the number 


of vertices, edges, and triangles of T respectively. The family {vi}i =denotes the vertices of 
the triangulation. 

For every k = 1,..., N t . we denote by (7}i=i,2,3 the vertices of the triangle T k eT, { vf, > 




the center of the edge linking v!f to Vj, and vfi = ^ J2i=l v i the center of mass of the triangle (see 
Figure [b]). Analogously we denote by {/f}o<i<3 the values of the function / at location {7}o^*^3 
of T fc . 
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Figure 6: Labels of various points in the triangle T*.. 

5.2 P 0 and Pi triangular finite elements 

Let us start with the following definition, 

Definition 5.1. For a given triangulation T, we denote by Po(T) (resp. Pi(T)), the set of functions 
that are constant on the interior of each triangle and null on their sides (resp. the set of the continuous 
functions that are affine on each triangle). 

The elements of Po(T) (resp. Pi(T)) are completely described by their values {f^)k at the center 
of mass (uq )k (resp. their values ({/,■'"}li$. 3 )fc at vertices ({ v i}l<K3)k) of the triangulation. Note 
that for / £ Pi(T), and each triangle T*, £ T 

3 3 3 

v/ G Pi(T), /o fc = M) = f(±v?/3) = E/(«?)/ 3 = E/*/3. 

2=1 2=1 2=1 

On the other hand, if / £ Po(7~), as / is null on the sides of each triangle, we can not calculate /g 
by the values at the vertices, and we can only set 

V / £ P 0 (T), f% = M). 

We denote by po the L 2 -projection of Pi(T) on Po(T) where for / £ Pi(T), the function po(f) is 
the unique element of Po(T) such that 

Po{f)(vo) = f{vo), (5.1) 

i.e. the operator po replace the affine approximation of a signal on each triangle by a constant 
approximation using the value at the center of mass. 

A basis for Pi(T) is given by the family ,-,n v with <pi £ Pi(T) and = 5ij (Kro- 

necker’s delta), for every i = 1, • • • , N v . Then, every / £ Pi(T) can be written as 

N v 

VxeT f(x) = fjVj( x ) i fj = f( v j) ■ 
j= i 

Remark that there exists a bijection between Pi(T) and M. Nv , defined by the following operator 

N v 

Pi = (/i, • • -Jn v ) ->• f = Yf j <P j Z Pi(T) (5.2) 

3 =1 

5.3 Discrete operators 

For every k = 1, • • • , iVj the area of the triangle T) ; is denoted by |T^| and is equal to ^ 11 riT k \ \ 
where riT k = (v\ — v\) A (v$ — v^). 
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5.3.1 Discrete norms on signal 

Depending on how the continuous signal is discretized, various methods may be used to compute 
the norm of the discrete signal. 

Formula exact on Po-finite elements Let p ^ 1 and f : T ->Ka function in L P (T, M). The 
p-power of the discrete L p -norm of / is simply defined as 

N t 

L p 0 [f,T] = J2\Tk\\fo\ P - ( 5 - 3 ) 

k =l 

Formula ( |5.3[ ) is exact for signals that are (almost everywhere) constant on each triangle, i.e. we 
have L%[f,T] = \\f\\ P L r {T) for any / G P 0 (T). If / is in Pi(T), then we have Lg[/,7] = \\po(f) \\ P LP(;r y 

Formula exact on Pi finite elements In work, we need to compute L p norms with p = 1,2. 
Thence, if we consider a signal / G Pi (T) , it is sufficient to use a formula exact on piecewise 
quadratic polynomial to fit our needs. For instance the Newton-Cotes formula for triangular mesh 
(see |2] pages 178-179) is 


L\[f,T\ = 


1 Nt 
6 k =i 


Tk |(|/i fc 2 


p +i/i fc 3i p +i/2 fe 3 r 


(5.4) 


where for a piecewise linear signal / G Pi(T) we have //) = |(/f + f p ) for any k = 1,... ,N t and 
1 ^ i < j < 3. 


The formula (5.4) is exact if p = 2 (i.e. L\[f,T\ = ||/|| 2 L 2 (j-\) since the function / 2 p is a 
polynomial of degree 2 for any k = l,...,Nt. The formula ( |5.4[ ) is also exact if p = 1 and / 
has constant sign on each triangle (i.e. we have L\[f,T\ = ||/||z,i(T) if / ^ 0 or / ^ 0). If the 
signal / G Pi(T) has a changing sign on a given triangle Tk, we briefly explain how the quantity 
L\[f lT , {Tfc}] is computed in practice. First, remark that at least two of the /f ’s have to be of the 
same sign. Without lost of generality, assume that < 0 and /| > 0, ^ 0 (as in Figure 0. As/ 

is linear on Tk, we can easily compute the coordinates of the points v\,v§ belonging to the edges of 
Tk and satisfying 

f(tv 4 + (1 — t)wg) = 0, for t G [0,1]. 

The triangle T/. can be split in three subtriangles T^, T k 2 , T k 3 with vertices (v\,v\,v§), (v 2 ,v^,v\), 
(v 2 ,V 3 ,V;() and with corresponding signals Pi (/i, 0,0), Pi(/ 2 , 0 , 0), P\ (/f, /a , 0) respectively (the 
bijection P\ : M 3 {Tk} is defined by (5.2)). Thence, the exact L 1 norm of f\ T G Pi({Tfc}) may 


now be computed by applying the Newton-Cotes formula on this three subtriangles: 


L\lfiT k , UU] = 3 E ir fc |(l/3l + l/SI + \&- 

1=1 

The general case follows as L\[f,T] = J2kh L\[f\ T ,{Tk}]- 


5.3.2 Discrete norm on signal derivative 

Let assume that / G Pj (T). We give in the following the discrete operators needed to compute 
the discrete H 1 and BV norms. For every Tk G T and for every / G Pi(T), the gradient of / on 7). 
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I/: 


v\ 


Figure 7: The L 1 norm on a triangle T k of a piecewise linear signal (/f, fif, ) is the non-oriented 
volume between the blue and green triangles. This figure illustrates how T k is divided in 3 subtri¬ 
angles to compute of the L 1 norm of a signal. 


can be computed by 


[Vr f]r k = 


e§ A e? 


fc "e§Ae §|| 2 


A (/M + / 2 fc e| + /|e|) 


(5.5) 


where 

In this framework, gradient V 7 -/ is constant on each triangle and by convention is null on the edges. 


e l — ^3 - v 2 > e 2 — V 1 ~ v 3 > e 3 — v 2 ~ V 1 


Total variation. The total variation of / 6 Pi on T is given by 


A7 


n/,n = Ei^ill[vr/]Tj. 

fe=l 


(5.6) 


7/ '-norm. For / E Pi(T), the square of the L 2 -norm of the gradient is given by 

N t 

H[f,T] = E l^~M||[Vr f]T k \\ 2 ■ (5.7) 

k= 1 

Remark that H[f,T] = To[||V-r/||,T] as / E Pi and by formula 

5.3.3 The fvarifold norm 

The fvarifold associated with the discrete fshape (7~,/) is naturally defined by However, 

in order to simplify the computation we approximate such a fvarifold by a discrete measure 
Such an approximation is set in the same way in the Po and Pi-finite elements setting. 

The discrete fvarivold norm is defined as 


(5.5) we have ||Vr/|| G Pq(T). 


N t 


v(T,f) - E man,® 

k =1 


(5.8) 


v k o = l(vt + v k 2 + v k 3 ), / 0 fc = f(v8), V k = Spanjuf - v*, v k 3 - . 


where 
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The discrete fvarifold-norm is defined by 

Var[/i(r,/)] = \\H( T ,f)\\w' 


(5.9) 


which can be easily computed by ( |2.3[ ). We point out that by formula (5.1) we have 

^(T,/) = ^(T, P0 (/)) V / G Pi(T) . 


Finally, as we identify the tangent space with the normal one, accordingly to (2.1), we define the 
discrete RKHS by the following kernels 

\y- v '\ 2 2 (l-<ny r V ,) 2 ) p-f'l 2 

k e (v,v') = e , k t (V,V') = e ,k f (s,s') = e , 

where ny denotes a unit normal vector to V. 


6 Discrete problems and T-convergence results 


We now come back to the continuous “functional” matching problem defined by formula (3.1). 


The aim of this section is to show the T-convergence of well chosen discretized versions of these 
problems to their continuous counterpart. To this end, we have to compare functionals defined 
on different spaces. This can be done by projecting the discrete functional on the limit surface as 
described in Definition 4.11 Following [14], this leads to define a suitable topology (called in the 


following S-topology) in order to set a definition of T-convergence. 

6.1 Discretization procedure 

First of all, we discuss how to discretize (on an admissible triangulation) a signal defined on a 
surface. Let X be a surface verifying Hypothesis [l] and T be a /i-admissible triangulation for X for 
some h > 0. We denote by / a signal on X belonging to L 2 (X) (resp. H 1 (X) or BV(X)). The aim 
of this Section is to define from / a Pq or P\ finite element /), on T. 


According to property ( ii ) of Definition 4.4 by using the map Ext defined in (9.13), we can extend 
the manifold A to a larger suitable manifold X +Tl such that T C A4(A' +,? ). Moreover, because of 
the signal / can be extended to a signal / defined on X +r> which is L 2 (X +V ) (resp. 


9.5 


Theorem 

(X +ri ) and BV(X +V )) and with a small norm on X +v \ X. 

Let T be a triangle belonging to A 7 -. Depending on the regularity of / we consider two dis¬ 
cretization procedure: 

• the Pi-element: the piecewise linear fh corresponding to / on T is defined as 

fh = Pl(f{K x +v{vi)), f{Tr x +v{v2))J{TT X +v(v3))) 


where vi,V 2 ,V 3 denote the three vertices of T and P\ is defined by (5.2). 


the Py-element: as every Po _e l emen t is nu ll on each side of T, we define the piecewise (ae) 
constant fh corresponding to / on T as 


fh = 



on the sides of T , 

vi)) + f(^x+v(v: 2 )) + f{n X +ri(v 3 ))) otherwise. 


where Vi,V 2 ,vs denote the three vertices of T. 
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6.2 Preliminary result 

In this section we prove a technical result that is useful in the following. 

Lemma 6.1. Let X be a surface verifying Hypothesis [l] and T a /i-admissible triangulation for X 
verifying Hypothesis [2] Then we have 

sup II P(T,f) ~ P(XJ*)\\w' = °( h ) , (6- 1 ) 

/eP 0 (T) 


and 

v/ € Pi(T) ||/* (r>/) - n {x ,fi)\\w = 0(h)( 1 + ||V/|| L 1 (T) ). (6.2) 


Proof. Let / E Po(T). By the change of variables x = irx{y), because of Hypothesis [2] and formula 
(4.5), we have uniformly in / E Pq(T) and ip E Cq(M 3 X G(3, 2) X M), 


lx 


<p(x,T x X, f(x))dU 2 (x) - / <p(y, T y T h \ f(y))dH 2 (y) 


Moreover, because of property (iv) 
/z. —>• 0. Then 


J T U 


in Definition 


£ 


4.4 


( MM(T,/) - P(x,f*)) 




L°°(R 3 xG(3,2)xR) 0(h ). 
Hff'outj-) converges towards the null varifold as 

L°° (R 3 x G(3,2) xM.)0(h) . 


Now, if / E P 0 (T), then / = /q on the interior of every triangle T k . Then, for every (p E 
Cq(M 3 x (7(3,2) x M), we get 

N t . 

\(P{T,f) ~ ^(T,/) )(</>) I ^ H / \<P( x ^ T xT,f(x)) - <p(v%,V k j£)\dH 2 (x) 

k= l JTk 
N t 

^ ^2 I|¥’IIc 1 (R 3 xG(3,2)xR) / ||(®, V k , f(x)) - (vq, Vjc, /o )||r 3 xG(3,2)xR d'h 2 (x) 

k= 1 JTk 

N t 

^ l^fclll < / J llG 1 (K 3 xG(3,2)xR)diam(Tfc) ^ , H 2 (Th) Hv :> llG 1 (R 3 xG(3,2)xR)C ) (/i) . 
k =1 

where IVt denotes the number of triangles belonging to T. So, by taking the supremum on the unit 
ball of Cq(M 3 x (7(3,2) x M) and by using (2.2), we get 

II d(TJ) - P(TJ)\\w C\\mrj) - /^(T,/) II(Cq(R 3 xG(3,2)xR))' ^ 0(h)PL 2 (T). 


Then, we have \\[ji(x,f e ) ~ P(TJ)\\w ^ 0(h) which proves ( |6.1[ ). 

If / E Pi (T) the proof is similar, but the bound of | (^(T,/) — ^(T,/))(^)l depends on / because 

Vx E T k | f(x) - / 0 fc | 11V/11 L i (Tfe )diam(Tfc) 

so that, by the same arguments, we get 

II P(T,f) ~ P(X,p)\\ W' ^ 0(h)(l + ||V/|| L i(7-)) . 


which proves ( 6 . 2 ) 


□ 
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6.3 The L 2 -model 

We prove a T-convergence result for the discretization of the problem ( |3.2[ ). The proof is given 
in the relaxed setting of measures defined in Definition |3.3| We refer to Section |3.1| for the main 
notations. Let {Th} be a sequence of /i-admissible triangulations for X verifying Hypothesis |2j The 
discrete energy E\ : Po (Th) -^KU {+ 00 } is defined by 

E h (h) = fiSlA, T h \ + ^Var[„ (7L/ ,j - (6.3) 

where 33 % is a /i-admissible triangulation for Y verifying Hypothesis [ 2 ] and < 37 ,, G IPo(34) is the discrete 
signal on 34 obtained by applying the discretization procedure described in Section |6.1| to g. We 
have in particular 

\\P(Y,g) ~ A t (y h , flh )llw’' —• > 0 as h — > 0 
We define also (for every h ) the following discrete problem 


inf 

fh&’ofjh) 


E h (f h ). 


(6.4) 


We also suppose that 7 // 7 W is independant from h and verifies the hypothesis of Theorem 3.2 for 


every h. An existence result for the discrete problem can be proved easily in the discrete setting. 


Proposition 6.2. For every h > 0, there exists at least a solution to problem (6.4). 


Proof. Simply notice that Eh is continuous in fh, that Lq[//i,7/(] = J2k =1 l^fellAC^o )! 2 i s a coercive 
bilinear form on P 0 {%) and that Var [p(T h , fh ) ~ P(y h ,g h ) 1 = II P(T h ,f h ) ~ P(y h ,g h )\\w' is bounded. 

□ 


As the discrete and continuous functionals are not defined on the same space, we consider first 
the S'-topology to compare the signals. We then define the T-convergence (see e.g. M) for the 
functional: 

Definition 6.3 (5-topology and T-convergence). Let A be a surface satisfying Hypothesis [IJ 
Let (Th)h be a sequence of /i-admissible triangulations for X verifying Hypothesis [ 2 ] and / a function 
defined on X. In this definition, (fh)h denotes a sequence of functions such that fh G Po(7h) for 
every h > 0 . 


Cf 

We say that (fh)h converges to / with respect to the 5-topology (fh f) if and only if 


P(xjf) A P(xj) in M x , 


(6.5) 


4.11). This defines a 


where for any h > 0, /( is the projection of fh onto X (see Definition_ _ 

generalized convergence between functions with different supports which is used to defines the 
T-convergence. 

We say that (Eh)h T-converges to E if: 

( i ) Lower bound: for every / G L 2 (X ) and for every sequence (fh)h such that fh /, we 
have 

E(f) ^ liminf E h (f h ) ; 

h —>0 
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(ii) Upper bound: for every / £ L 2 (X ) there exists a sequence {fh)h such that fh —^ /, and 


U(/) ^ limsupU/ l (//i). 
h- 5-0 


Theorem 6.4. The sequence (Eh)h T-converges to the functional E defined in (3.2) with respect 
to the S'-topology. 

Proof. Lower bound. Let (fh)h be a sequence of functions such that fh £ Po(7h) for every h and 
fh f £ L 2 (X). As the functional E defined in (3.4) is lower-semicontinuous with respect to the 
weak convergence of measures, we obtain 


E U) = E iV(x,f)) < liminf E(ti(xjt)) 
Now, by change of variables and because of (|2]) and (|4.5[) we get 


L\ 


o ihM = \\f h \\i HTh) > mhtfr) > ( i+om\fh\\h {x )• 


Moreover, from Lemma 6.1, we have 

Varner,,,A) ~ V(y h ,g h )\ = II P(T h ,f h ) ~ P(y h ,g h )\\w' = \\E(x,fD - P{Y,g e h )\\w' + 0(h) 
Finally, by ( |6Ul ) and ( [(k8| ) and since || H(y, g ) ~ H(y h . S]h )\\w' = o(l) 

Eh(fh) ^ (1 + 0(h))E(fff) + o(l) ^ (1 + 0{h))E{n^ X ftf) + o(l) 


( 6 . 6 ) 

(6.7) 

( 6 . 8 ) 

(6.9) 


and, because of (6.6), this proves the lower bound condition. 

Upper bound. As C' 1 (A) is dense in L 2 (X ) we can suppose that / is smooth (see Remark 9.6). 
We can then conclude by a diagonal extraction argument (see below). 

For every sequence ( Th)h of /i-admissible triangulations for X verifying Hypothe sis [2] we can define 
the functions fh £ Po (Th) defined on Th by the Discretization Procedure of Section 6.1 for the setting 
of the Fo-fiaite elements. Then, we define a sequence fh. £ L°°(Th) such that sup^ \\fh\\L°°(T h ) ^ 

Definition 


L°°(x)- Thus, because of ( iv) 


4.4 


we obtain 


\\fh\\L 2 (T° ut ) — 0{h) 


Now, we get f^—^f strongly in L 2 (X), and, because of (5.3), Proposition 
[2] by change of variables, we obtain 


4.10 


( 6 . 10 ) 
and Hypothesis 


WfhWl^x) = L oifh^ h } + 0{h) 


so that 


Lo[fh,Th] —> ||/Hl 2 (x) as 0 . 


In particular, (up to a subseque nce) f[ -4 / a.e. which implies that fh — 1 / (meaning that 


H( X —*■ //(x,/)) and, since Lemma 


6.1 


and Lemma 


3.1 


we also get 
2 


Va x \P{T h ,fh) ~ P(y h ,9h )] - ll^(Th,/h) _ P(y h ,g h )\\w' -4- IIM(X,/) - P(Y,g)\\w' 


Then, we have proved that fh —*■ /, and 


E(f) = lim E h (f h ). 
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Diagonal argument. Let (/n)neN be a sequence in C' 1 (X) converging to / G L 2 (X) in L 2 
and pointwise. As the fvarifold-norm is continuous with respect to the pointwise convergence of the 
signal we have 

lim E(f n ) = E(f). 

n—>-+oo 

Previous proof shows that, for every n G N we can define a Po~ e l emen ts sequence p n (formerly called 
fh) such that 

VneN, .lim E(f n ) = E(f n ). 

2—>-+00 

This means that 


Ve > 0, 3N(e) Gff:n)]V4 | E(f n ) - E(f)\ < e /2 


and 


Vn € N, 3J(n) £N :i) /(n) =7 | E(f n ) - £(/ n )| < l/n. 

Let e > 0, the sequence (/n^) n eN verifies 

n ^ iV(e) =► | E(f^) - £?(/)| ^ | E(fW) - E{f n )\ + \E{f) - E{f n )\ < l/n + e/2 
which yields to 

Vn ^ max{7V(e),2/e} =7 \E(f^) - E(f)\ ^ e. 

The last inequality is valid for arbitrarily small e > 0 and it eventually proves that 

lim E(f^)=E(f). 

n—>•+oo 


m 


For the following result, we consider the assumptions of Theorem 3.2 that is that A is a C p 
oriented connected compact 2 dimensional submanifold (p ^ 2 ), Y is a finite volume 2 -rectifiable 
subset of M 3 and that W is continuously embedded in Cp{M. 3 x G(3,2) x M) with m ^ p. 


Theorem 6.5 (Convergence of minima). Assume that hypothesis of Theorem 3.2 hold and 


that X satisfies Hypothesis [T| Let ( Th)h be a sequence of admissible triangulations for X verifying 
Hypothesis [2j Then if 7 / p/w is large enough, 

lim min Eh(f)= min E(f). 

h^O /SPo(A) f£L*(X) 

Moreover, if (//) is a sequence of minimizers of the functionals Eh, then, up to the extraction of a 
subsequence, there exists /* G C' P_ 1 (A) such that fh k 

Proof. The inequality 


f* and £(/*) = min /ei 2 ( X ) E{f). 


lim min Eh(f ) Sj min E(f). 
h^o /ePo (T h ) f&L2(X) 


( 6 . 11 ) 


follows easily from the upper bound condition of T-convergence. 

Let us look at the other inequality and consider fh G PoiTh) be a minima of Eh for any h > 0. 


Because of Theorem 3.2 (in particular (3.3)) we have 


\\f h \\ L ~ {Th) ^c™{u\T h ) + u\y h )) 


if 
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As the areas of the triangulations converge to the areas of the respective surfaces, we can assume that 
the sequence {f^)c is uniformly bounded on X. It yields that the measures \i^ x ^ are supported on 
a compact set. Thence, the sequence {g^ x f e ) } i s tight and, because of the Prokhorov’s theorem, it 


weak-* converges (up to a subsequence) to some £ M x . Then, by using the lower semicontinuity 
of E with respect to the weak-* convergence of measures and the fact that E is minimized by a L 2 - 
function (see Section 3.1), we get 


, -ECf) ^ ^(Moo) < iiminf E(g {xft) ). 
feL 2 (X) h-± 0 v 


Now, since (3.5) we have E(g^ x ^) = E(ffc) and in the same way we get (6.9) we have (1 + 


0{h))E(fl) + o(l) 4 E h (f h ). It yields, 


min E(f ) ^ E(u 00 ) ^ Iiminf E(u (x ^ Iiminf min Ehif ) 

f&L 2(.Y) ' h^o ^ h^O / G Po(T h ) 


that gives the needed inequality. Using (6.11), we get also that minj eL 2 m E(f) = E(fi qq), and 
from Theorem 3.2 that there exists /* G C P ^ 1 {X) such that = g(x,f«)- Thus, fh /*• 

□ 


6.4 The 4>U-model 


We study in this section the T-convergence of the discretization of (3.7). The Sobolev case can 


be treated by the same arguments by considering the discretized H 1 -norm instead of the £>U-one 


(see Remark 6.10). In this section we suppose that {Th} (resp. {34}) is a sequence of /i-admissible 
triangulations of X (resp. Y) verifying Hypothesis [ 2 ] and belonging to Tc for some C > 0 (see 
Definition 4.2). Then, for every h, we define the dicrete energy E^ : Pi (7^) ->IU {+ 00 } by 


E h (h) = (.+ V[A,7i])) + ^Var [ HTh:Ih) - 


( 6 . 12 ) 


where G Pi (34) is the discrete signal on 34 obtained by applying the discretization procedure 


described in Section 6.1 to g. That defines (for every h > 0) the following discrete problem 


inf E h {f h ). 

fh&i(Jh) 


(6.13) 


Proposition 6.6. For every h > 0, there exists at least a solution to problem ( |6.13[ ) . 

Proof. For every / G Pi(7fc) we have ||Vr h /||Li(r h ) = v [f,%] and = \\f\W(T h )- Then 

every minimizing sequence is bounded in BV so that it weak-* converges in BV and a.e. (up to a 
subsequence). The result follows from Lemma |3.1| the lower semicontinuity of the total variation 
with respect to the weak-* convergence. 

□ 


As in Section [6.3[ we introduce the S-topology needed to define the T-convergence for functionals: 

Definition 6.7. [S’-topology and T-convergence] Let X be a surface satisfying Hypothesis [I] and 
C > 0. Let (' Th)h be a sequence in T c of admissible triangulations for X verifying Hypothesis [2J and 
/ be a function on X. In this definition, (fh)h denotes a sequence of functions such that fh G Pi (Th) 
for every h > 0. 
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We say that (fh)h converges to / with respect to the 5-topology (fh —^ /) if and only if 

I™ Wfh ~ /IIli(a-) ( 6 - 14 ) 


where for any h > 0, /)( is the projection of fh onto X (see Definition 4.11). 

We say that ( Eh)h T-converges to E if: 

(i) Lower bound: for every / E BV(X) and for every sequence (fh)h such that fh —*■ f, we 
have 

£(/) ^ liminf E h (f h ) 
h—>0 

(ii) Upper bound: for every / E BV(X) there exists a sequence (fh)h such that fh ^ f 
and 

E(f) ^ lim supEh(fh) 

h->0 


Theorem 6.8. The sequence ( Eff)h defined by formula (6.12) T-converges to E with respect to the 
5-topology of Definition |6.7[ 


Proof. Lower bound. Let (fh)h be a sequence of functions such that fh E Pi(7h) for every h > 0 
and fh / E BV(X). Thus we have 


where /)( is the projection of fh onto X (see Definition 4.11). Without loss of generality we can 
suppose that 

sup E h {fh) < oo . (6.15) 

h> 0 

By Proposition 4.12 Hypothesis [2j and condition (iv) in Definition 4.4 we get 

\\fh\\BV(T h ) ^ (1 + 0(h))\\fh\\ B v(X), 

which implies 


(l + om\fi\\ BV{X ) < \\fh\\ B v{T h ) = L \[fh, Th\ + V[fh, Th\. 

Moreover, by Lemma 6.1 and equation (6.15), we have 


W^iThJh) ~ l*{x,f*)\\w' < °( h ) 


(6.16) 

(6.17) 

(6.18) 


As /)( f. we get that 


Vai '[^(T h ,/ h ) - A*(y h ,sfc)] ->• ll^(X/) “ HY,g)\\w' ■ ( 6 - 19 ) 

Now, the HU-norm is lower semicontinuous with respect to the L 1 -convergence then, because of 


(6.17) and (6.19), we get 


7 w 


E (f ) < bminf ( Wfh\\BV{x) + —Var[/i (7 - hi/h) - A*(y h , Sh) ] ) < HminfU ft (/ h ) 


h—>0 


( 6 . 20 ) 


Upper bound. As every function in BV(X) can be approximated by a sequence of C 1 (X) 
functions with respect to the strict topology (see Theorem 9.5), we can assume that / E C 1 (X). 
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The general case ensues by a diagonal extraction argument similar to the one given in the proof of 
Theorem 16.4[ 


Section 6.1 we detailed the way we can define fh £ Pi(7ft) for every /r-admissible triangulations 
7ft for X. We remind that, in the discretization procedure, for every h > 0, we define an extension 
X +ri of the manifold X and, afterwards, we extend / to a function / £ BV(X +ri ) that is used 
to define fh- Of course the extension X +ri can depend on h, but in the following we prove some 
estimates for a fixed h > 0 and only write X +v . In particular, we have 


W 1 ’°°(X+ r i) 


€ 




( 6 . 21 ) 


which implies that a fh € W 1,00 (Th) and swp h ||/ft||jyi,°o( 7 ', i ) ^ WfWw 1 - 00 ^)- This is essentially due 


to the extension procedure by symnretrization defined in Theorem 9.5 


Since Theorem 3.1.6 p. 124 in [9] (applied with k = s = 0, p = oo, m = q = 1), we get 

11/ O 7 T X +ri — fh\\ bv (Th) = 11/ ° ^x+n — fhllw^iTh) ^ Ch ' 2 \\f ° ^ X+Alw 1 ’™ (T h ) (6.22) 

where C > 0 is a constant and / is the extension of / on the extended manifold X +ri . More precisely, 
to obtain the previous estimate from Theorem 3.1.6 p. 124 in j9], we need to sum on all the triangles, 
use (4.2), and the fact that the area of every triangle is bounded by n(h/2) 2 . 

Now, because of ( iv ) Definition |4.4| we have 

ll//i|lsy(r° ut ) = O(h ), (6.23) 

so that, by Proposition |4.12~ and ( |6.21[ ), we have 

11 / “ /ftllw(X) = 11/° Tv +r ! - fh\\BV{T^) + 0(h) , 

11 /° TTx+vWw 1 ’™^) ^ WfWw^^x) + 0(h) 

that, with (6.22), give 

11 / ~ /ftllw(x) = 0(h). 

This proves in particular that /^ —>• / strongly in L 1 (A') (which means that fh f) and the 
convergence of the total variation term, because 

WVxfhWmx) = II’ VtMlw?) + 0(h) = V[f h , Th] + O(h ). 

Moreover, Proposition |4.12| implies 

ll/ftlUpx) = WfhW^iTh) = L\[fh,T h \ + 0(h), 

so that 

L\ [A, 7ft] —>• ||/||^i(^) as h -> 0. 


□ 


The convergence of the varifold term follows by Lemmas | 6 .1| and 3.1 Then 

E(f) = lim Eh(fh)- 

h—>0 


Theorem 6.9 (Convergence of minima). Let (7ft)ft be a sequence of admissible triangula¬ 
tions for X verifying Hypothesis | 2 j Let (fh)h be a sequence of minima of Eh (he., Eh(fh) = 
min^gp^T^ Eh(f)). Then, (fffjh weakly-* converges in BV(X) (up to a subsequence) to a minimum 
of E and 

lim min Euif) = min E(f). 
fc->o/ePi(T h ) feBV(X) 
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Proof. From the previous proposition we have that Eh T-converges to E. We consider the sequence 
(fh £ Pi(7ft)) of the minima of Eh and, without loss of generality, we can also suppose that 

sup E h (f h ) < oo . 

h> o 


As for (6.17), we have 


II f e h\\BV(x) < L\[f h ,T h ] + V[f h ,T h \ + 0(h), 

so that {fh)h is uniformly bounded in BV(X). Then, there exists a su bse quen ce (not relabeled) 
such that ( fh)h weak-* converges to f°° in BV(X) and, by Lemmas 6.1 and 3.1, we get 


Varft[M(77,A) ~ P(y h ,g h )\ \\P{xj°°) - V(Y, g )\\w'■ 


Then, similarly to (6.20), by lower semicontinuity, we get 


min E(f) ^ E(f°°) ^ lim min Eh{f). 

fesv(x) h^ofeMTh) 

The other inequality follows from the upper bound condition of T-convergence. It follows in particular 
that f°° is a minimum of E. □ 

Remark 6.10. (The Sobolev case) The discrete problem for the H 1 -model is defined in the same 
way by considering the discrete energy 

(£?[/, 71 ] + + ^Var k [n nJ) - 


Then, Theorem |6.8| and Proposition |6.9| are still true in this case and their proofs can be straight¬ 
forward generalized by considering as S'-topology the L 2 -convergence of the projection (i.e., fh~^f 
if fl f strongly in L 2 (X )) and using the compactness with respect to the weak topology of H 1 


(see Remark 3.6) 


7 Numerical results 


In this section we show some numerical results illustrating the effect of the different penalty 
terms based on L 2 , H 1 or BV norms on the optimal solution. All the codes and data used in this 
examples are available at https://github.com/fshapes/fshapesTk 

The energy is discretized in the framework of the finite elements as detailed in Section[5]and 6.1 


Note, that all the simulations are performed with a number P of points fixed. Moreover, in order to 
apply a gradient descent algorithm we choose to consider a smoothed version of discrete BV norm 
(the discrete L 2 and PL 2 being differentiable). For the numeric simulations, we define the discrete 
operators (5.6) with the smoothed norm ||-|| e = \J (-) 2 + e 2 instead of the classic euclidean norm. 


This makes differentiable the L 1 -norms and the evaluation of the jump part of the total variation, 
for every finite element. Note that such a smoothed energy T-converges to the initial one as e goes 
to 0, and may justify such an approximation. 











(a) Source (X, 0) 


(b) Target (Y,g) 


(c) Relative position (mesh view) 


Figure 8: 8a the initial fshape is a fat smooth star (same as Figure la). 8b: the target fshape (Y. g ) 


is a thin star with Lena picture as grey level signal (black is 0 and white is 1). Compared to Figure 
lb the pixel effect on the target has disappeared because of the Pi discretization (see Section [5]). 
The blue cylinder is just drawn to help to see the perspective. 

7.1 The melting smooth stars 


We consider the matching problem between the two fshapes drawn in Figures 8a and |8b[ These 
fshapes are the Pi versions of the fshapes pictured Figure [lj Recall that we are interested in finding 
a signal / defined on X so that the fshape (X, /) is as close as possible to the t arget fshape ( Y,g ). 
The geometric support X remains fixed. As the P 2 -model is given by Example 3.4 we show here 
only the effects of the BV and H 1 penalties. 

We analyze in particular the impact of the volume and gradient term by considering the following 
weighted energies of the residual signal: 


E(f) =pen(f) + ||| m x ,f) ~ V(Y,g)\\w'’ 


with 


1 1 \ _ ) “11/llppO + P\Df\{X), 
peni/, '\|b"2 ' «"V7..rl|2 


BV model 
H 1 model 


ILRX) + P\\V xf\\ L 2( X )i 

For sake of readability we have written the continuous version of the functional E but the mini¬ 
mization is obviously performed on the discrete version. The optimal fshape is defined by ( X,f *) 
where f* is a (possibly local) minimum point (of the discrete version) of E found after an adaptive 
gradient descent. 

In Figure [9] and 10 we show some simulations for different values of the parameters a > 0 and 
f3 > 0. The main difference with Figure [2] is the absence of oscillations for the signal f* in the non 
overlapping part of the source and target fshapes. Nevertheless, if /3 is chosen too small, oscillations 


may appeared as in Figure 10c 


7.2 The OCT dataset 

We consider the two fshapes shown in Figure These data are obtained by segmentation 
from volumetric images of the retina acquired by optical coherence tomography ( |T6| ITT]). The 
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(a) The target (Y, g) (b) a = 0.00125 and 0 = 0.125. 



(c) a = 0.00125 and 0 = 0.0125. (d) a = 0.125 and 0 = 0.0125. 


Figure 9: Results for the BV penalty: en face view of the target (Figure 9a) and of the optimal 
signals f* defined on X (Figures 


9b 


9c 


and 


9d). 
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(c) a = 0.001 and 0 = 0.01. (d) a = 0.001 and 0 = 0.001. 


Figure 10: Results for the H 1 penalty: en face view of the target (Figure 10a) and of the optimal 
signals f* defined on X (Figures 


10b 


10c 


and 


lOd). 
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surfaces represent the lower boundary of the Nerves Fiber Layer (NFL) and the signals represent 
the thickness of the NFL. The variation of the thickness of the NFL may be an indicator of glaucoma. 
Considering Figure 11c or 12a (same scene from two different points of view), one can clearly see the 


misalignment of the two fshapes. This problem is very common in image analysis and the example 
of the smooth melting stars was a way to illustrate this problem. 




(a) Source data (X,f). 


(b) Target data (Y,g). (c) Relative position. 



Figure 11: En face view of 2 particularly non-align individuals of the OCT dataset. 


In practice, the usual transformations acting on fshapes are both geometric and functional as 
described in [Bj. We then show some numerical results for the energy 


E{f,p)=pen{f) + pen(p) + A(p ■ (X, f 0 + f),(Y, g)), 


where 

• pen(p) is the energy of the deformation p : M 3 — y M 3 sending the source fshape on the target 
fshape. We use the Large Deformation Diffeomorphic Metric Mapping framework for fshapes 
(see 0 ) to generate such a deformation and to compute its energy. 

i a \\f\\ 2 L\x) L 2 model 

• pen(f) = | a||/|| L i (x) + P\Df\(X) BV model. 

[ a \\f\\h{x) + P\\ v xf\\h( X ) Rl model 

The signal / : X —>• M was called “functional residual” in [6] and is added to the original signal 
/o defined on X. The parameters a and (5 are positive constants. 

• the attachment term is 

■ (X , /o + /), (Y,g)) = 7 w (\\P(v(X)j 0 +f) ~ H{Y, g )\\w' + ll/^MdATo) ~ ^(dY^Ww 1 ) ■ 

It is a sum of two fvarifold terms: the first one involves the whole fshapes and the second 
one involves only the boundaries of the fshapes. This latter term enforces the matching of 
the boundaries and no signal is considered on the boundary. The parameter yw is a positive 
constant. 


The optimum value for the deformation and the signal found after a gradient descent on E are 
denoted by p* and f* respectively. The optimization in p and / is performed simultaneously in 





(a) t = 0 


(b) t = 1/5 


(c) t = 2/5 



(d) t = 3/5 (e) t = 4/5 (f) t = 1 


Figure 12: Perspective view of the geometric-functional shooting of the two fshapes presented Figure 

El 


order to find a deformation ip* taking the signal into account ( i.e. matching loci with similar signal 
patterns together). In Figure 12 we show the results of the numerical experiment with the H 1 
penalty. 


In Figure [T3j we show some results for the different models with a common 7 w and the best 
parameters a, /3 we were able to find for each model. We represent only the signal f* found after 
the gradient descent and we do not discuss here the quality of the geometric matching. We can 
observe oscillating signal with the L 2 model. Regularizing the gradient of the penalty on / improve 
the quality of the result. Nevertheless, some low frequency oscillations are still present on the non 
overlapping part. The results with BV and H 1 model are much better as no oscillation may be 
observed. The H 1 model may be preferred on this example as the signal f* has higher values. 


8 Conclusion 


Previous examples show the effects of a BV or H 1 penalties on the matching between fshapes 
and point out the quality of the results compared with the L 2 model. This is due essentially to 
the behavior of Gaussian kernels which can make the energy lower if the signal is oscillating (i.e., 
|/i — / 2 I large). The BV or H 1 models prevent from the oscillations and avoids uninteresting optimal 
solutions. 

The paper takes into account the model with fixed geometry. A further problem concerns the 
study of the effects of the BV or H 1 penalties on the geometric evolutions. In the simultaneous 
optimization of the signal and geometry, the L 2 -model can be trapped in poor optimal geometric 
configurations. The effects of the BV or H 1 models on these geometric phenomena has to be studied. 
The T-convergence for the global (geometry and signal) model is also to prove. These questions will 
be addressed in a future work. 
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(a) L 2 -model 


(b) L 2 -model (with regularized gradient 
descent) 



(c) BV model. (d) H 1 model. 

Figure 13: Optimal functional residuals (X, /*) with the various penalty terms. 

9 Appendix: functions of bounded variation on manifolds 

9.1 Definitions and main properties 

Here we consider a slightly more general framework where X is not supposed to be a finite 2D 
compact submanifold on M 3 but a general compact d dimensional manifold. The definition of BV 
functions is however dependent of the introduction of a divergence operator (or equivalently of a 
volume form) and a local notion of length, two things that are provided by a Riemannian structure. 

So, let X denotes a oriented smooth (at least C 1 ) compact d-dimensional Riemannian manifold 
possibly with boundary denoted OX and let volx be the associated Riemannian volume form. The 
boundary dX is supposed to be a C 1 compact (d— l)-dimensional manifold and we have volx(dX) = 
0. Finally, let us denote 

X Q = X\dX 

which is a C 1 manifold without boundary. When X is without boundary, the previous construction 
gives Xo = X. We say that / G L 1 (X, M) is a function of bounded variation on X if 

\D x f\(X) = sup | /div. Y (w)vol.Y | u G xj Ao), IMloo ^ l| < oo 
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where Xc(^o) denotes the set of C 1 vector fields u : X —v TX on X compactly supported in Xq and 
divx is the divergence operator on X defined by 


d 

divx(u) = g( e i, du(ei)) 

i=1 

where (ei, • • • , e^) is an orthonormal frame on TX. Here ||z/||oo = su Piex dx{u(x),u(x))2 where g is 
the metric tensor associated with the Riemannian structure. We recall the integration by part given 
Stokes formula : 

/ /idivv(ri)volY = — / u{h)xolx = — / g(Vh,u)v oIy (9.1) 

Jx Jx Jx 

where h G C\(X o) and u(h) denotes the derivative of h along the vector fields u defined by [u(h)] (x) = 
d x h(u(x)) for any x G Xq. Obviously we retrieve our previous submanifold setting when we consider 
the metric induce on the submanifold by the ambient space M 3 . In the other way, a much larger 
setting has been developed in a recent paper by Ambrosio et al. [1] that goes beyond our framework. 
The functional space BV(X) endowed with the norm 

\\f\\BV(x) = \\f\\mx) + \Dxf\(X) 


is a Banach space. 

Definition 9.1. The space BV(X) can be also equipped with the following convergences, both 
weaker than the norm convergence: 

1. Weak-* topology. Let {/«} C BV{X) and / G BV(X). We say that the sequence (f n )n 
weakly-* converges in BV(X) to / if 


fn L ^? f and D x f n 


Dxf , as n —> oo; 


where —*■ denotes the weak convergence in the space of measure on X. 

2. Strict topology. Let {f n } C BV{X) and / G BV(X). We say that the sequence ( f n ) n 
strictly converges to / in BV(X) if 

fn L -^ ) / and \D x f n \{X) —■> \D x f\(X) , as n-> oo. 

Remark that the following distance 


d(f, 9 ) 


\\f-g\\mx)+ \D x f\(x)-\D x g\(x) 


is a distance in BV(X) inducing the strict convergence. 

We just recall that the strict convergence implies the weak-* convergence (converse is not true). 
For more properties of functions of bounded variation we refer to H3 E]. Functions of bounded 
variation on manifolds are studied in [3]. In the following we remind the main properties used in 
this paper and we adapt the proofs to the manifold case if necessary. 

Proposition 9.2. A sequence (f n )n in BV(X) weakly-* converges to / in BV{X) if and only if 
( fn)n is bounded in BV(X) and converges to / in L 1 (A). 


Proof. This is Proposition 3.13, page 125 in [3 . 


□ 
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Finally, we remind the main compactness result for functions of bounded variation. 

Theorem 9.3 (Compactness). Let ( f n ) n be a sequence of BV(X ) such that ||/n||w(x) is un i~ 
formly bounded. Then {f n } n is relatively compact in BV(X) with respect to the weakly-* conver¬ 
gence. 

Proof. The proof is similar to the classical-one (see Theorem 3.23 page 132 in [3]) and can be obtained 
by using a partition of the unity and the extension theorem like in the next section. □ 


9.2 Approximation by smooth C 1 (A 0 )-functions 

In this section we establish an approximation result of BV functions by smooth functions. Inter¬ 
estingly, even if there exists several density results of smooth function using smoothing through the 
heat kernel semi-group for geodesically complete Riemannian manifold (see j 1 8j , [5J and USD, to the 
best of our knowledge, an approximation result in the case of a compact manifold with boundary 
does not seem to be available. We give below such a result : 


Theorem 9.4. Let X be a orientable Riemannian compact manifold with boundary OX and let 
Xq = X \ dX. For any h G BV(X, M) there exist a sequence h n G C' 1 (Aq,M) such that 


K -4 h Uo in L 1 (A" 0 , M) 
and 


\D x h\(X) = lim f |VL n |volx ■ 

n^oo J Xq 


(9.2) 


Our proof shares similar ideas than the proof of the classical approximation result in the case of 
an open set 12 G (see Theorem 3.9 in $)■ However, the introduction of a non constant volume 
term when considering the problem in local charts introduces several new elements that need to be 
considered. 


Proof. Let {Uf)i^i be a finite atlas on Xq and for any 1 ^ i ^ n let ipi : Ui —>• Vi be a local chart such 
that Ui is compact and tpi is the restriction to Ui of a C 1 diffeomorphisms from Ui —> Vi (such an 
atlas exists since X is compact). Let be a partition of unity such that supp (r]j) is compact 

for any j ^ 0 and there exists a partition (Jj)i^ n of N for which supp(ryj) C Ui for any j G J* 
which is locally finite on any Ui (i.e., for any x G Ui , there exists an open set Ui(x ) C U t such that 
supp (rjj) FI Ui(x) = 0 for any j G J* except on a finite number of j’s). 

For any j G Jj we consider £j such that d(ipi(supp(jif)). Vf) > £j and for e = (£j)j^o we consider 
the linear operator L e : BV(X, R) —> C 1 (Ao,M) defined by 

n 

L eh = Y J Pi(Y, (HO * %) 

i= 1 j£Ji 

where -0j : Vi —>• Ui is the inverse mapping of Note that we are using here classical notation of 
differential geometry for pullbacks where for any function £ G C c (C/j,M), ip*£ = £ o ifi and for any 
v G x;i(t/j), r ip*v = (difi)~ l v o ifi. Eventually, on every Vi, we introduce a t d/x the pullback of volx \ b ,. 
by ifi on Vi such that for any £ G C c (Ui ), we have f v . l\o\x = f Vi (iff£)aidx. 

Let 5 > 0. We can assume that for any 1 0 ^ n and any j G Jj, we have £j small enough so 
that 

[ I(H) * Ptj - (■ hr)j)\oiidx < . 

J Vi 
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n 

Since {i)j)j^ o is a partition of unity, we have h = X) h-Vj anc i 

i=i ieJ., 


[ \L e h - h\vol x ^ J5 [ 55 IV’i (^j) * Pej ~ '(p*( hr lj)\ a i dx ^ 2( 5 • (9.3) 

Jx i =1 “W jeVj 

This first inequality is enough to prove an approximation result in a L 1 sense. We turn now to the 
control of the total variation part. 

Let u G Xc(Xq). We have the following decomposition using the integration by part formula (9.1) 
for equality (a) and the classical integration by part on for equality (6) 


L 


ix 


n . 

-/?div_Y(w)volx = L ME ipi (hijj) * p £j ^dW x {u)volx 
i= 1 JX j&Ji 

“ “ 55 / U [$ ( 51 (H’) * Pej ) ) vol A- 

i=l JX V j&Ji J 

n « 

= -55 /. 55($■«)(55 ^ (HO * p£j) a i dx 

i =i je Ji is Ji 

- 55 / 55 [V>*(H')l * /5 ej -div(oi^«)da: 

i=i j- e j, 

n „ 

= 55 /. 55 *pe j )dx 

„■_ 1 J V 5 .V T 


i=i ‘ /yi je Ji 


= 55 /. 55 v , *(M div (^j([«iV’» */ 0 e 3 -))rf® 


»=1 jeJj 


- 55 / 55 V’t* (M (Mi«] * Pej ) (V’i 7 h ) dx 


Bj 


with 


■A-ij — 


'Vj 


55 V>*(M div (atiipi rjj ( (ip*u) * p £ .)) + di v(^*Pj(M*«] * % “ OilM*") * PeJ)) 


jeJi 


dx. 


However, for any <5 > 0, denoting | • | x the norm at i £ I induced by the metric we have for 
x G supp(r/j) 


= (d a fpi) _1 ( / d v ,. (v3 . (a;) _ y) ^u(^(<Pi(x) - y))p ej {y) d y) 

JVi 

{dx^Pi} 1 ( (dMviW-y)^ ~ dxP^uiMwix) - y))Pej{y)dy 


^ 1 + 

^ 1 + d 
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for £j small enough uniformly in u such that Hujloo ^ 1 and j G Hence 

E / E Vf(/»)div(ai^7/j((^ t *«) * p ej .))ds ^ XI / E r h( 1 + <5)d|T> x /t| 


ieJi 

^ (l + ^l^/iKX). 


(9.4) 


Moreover, we can assume for Ej small enough we have for x G supp(r/j) 
[a i V , *«] * P £ , - ai[(^i“) * Pei 


<Pi- 


OLi 


d x ip. 


-l 


ai(<pi(x) - y) - ai(ipi(x)) 


I Vi 


ai(<Pi(x )) 


dij>(<p i (x)-y) t PM'<Pi( i Pi(x) ~ y))Pej(y)dy 




so that 


E / 2 ^div f a i^i Vj 

i =i 7r •/. V 


[aiV’l^] * Pei - OCi[(lpiU) * p £j 


OLi 


dx 


<E/ 


* ./ 


Thus we have 


E^' < |£>x/i|(X)(l + 2<5) 


2=1 

Let us consider now the B^s. We have 




H WV, 


’Vi 


(aiip*u)(x - y)p £j {y)dy , (V*^(-0*^)) (x)Jdx 
{oii^*u){x), j - y)p £j (y)dy\dx 


[ a i rM*u{r i y j )dx+ [ {(a^*u)(x), (V^V^)) * p £j - ^*hV(^ Vj ))dx 

JVi JVi 


B 1 


Concerning the Bjj terms, we have 


B 2 . 


(9.5) 


(9.6) 


(9.7) 


E E B h = E L E ^“(Pj)volx = j v hu( l)vol.Y = o. 

i=i jeJ,- 


i=l JX j&Jj 


lx 


(9.8) 


For the terms, let us notice that supy. |o:j'(/^u| < oo uniformly in u (since ||u||oo ^ 1) and 
(ip*hV(ip*r]j)) * p £j —> (ip*hV(ip*r]j)) in L l (M. d ,dx) so that for Ej sufficiently small, we can assume 
that | Bfj | ^ (52 _ ( J+1 ). Summing along the indices, we get 


EE4 

i=ijeJi 


^5 


(9.9) 


and with (9.6), (9.8) and (9.9) we get eventually that for sufficiently small values of the Ej’s we have 
uniformly in u £~)q(X q) satisfying Huljoo ^ 1 


lx 


L e hdivx(u)volx 


^ \D x h\(X)(l + 25 + 8 2 ) + S. 


(9.10) 
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Taking the supremum over such u, we get 

\D x L e h\{X) ^ \D x h\(X)(l + 28 + 8 2 ) + 8. (9.11) 

Since 8 is arbitrary we have shown that there exists a sequence (£k)k^o such that L £k h £ C'°°(Xo), 

limsup \D x L £k h\(X) < \D x h\(X) . (9.12) 

k 

Moreover, since we can assume that L £k h —> h in L l (X ), we get from the lower semi-continuity of 
the total variation in L 1 that L £k h —> h for the S'-topology. □ 


9.3 Approximation by C ,1 (A")-functions 

We can improve the previous result by showing an approximation result by C' 1 (A)-functions. To 
this end we define an extension of the manifold and we extend the function / to such a larger manifold. 
Then, we can obtain the needed approximation by a suitable approximation of the extended function. 
We give the details in the following. 

In the following X denotes a C p (p ^ 2) compact oriented 2 dimensional submanifold of M 3 
with non empty boundary denoted dX. We denote n Xo the C p_1 vectors field of positively oriented 
normal along Ao = X \ dX. Note that n Xo can be continuously extended on X and denoted in 
that case n x . Thence, for any x £ dX we can define the unit vector v{x) pointing outward and 


orthogonal to both n x (x) and dX. For sufficiently small r > 0, the subset J\[ t {Xq) defined by (4.3) 
is a open set on such that the mapping 


V* x 0 : (x,t) ^ x + tn xQ (x) 

is a C p_1 diffeomorphism from Aox] —r, r[ to J\f r (Xo) so that considering its inverse = (7 r Xo ,t Xo ) 
we have 

2 = irx 0 (z) + t Xo (z)n Xo (ir Xo (z)) 

for z £ N r (X 0 ) and t: Xq (z) can be interpreted as the orthonormal projection of z on Ao and |t.Y 0 (- 2; )| 
as the distance from z to Ao- Now there exists ?/o > 0 small enough such that the mapping 
Ext : dX x] — 770, ??o[x] — r, r[— >• M 3 such that 


Ext (a;, s, t ) 


xo + tn Xo (xo) with xo = vrx 0 (x + sv(x)) if s < 0 
x + su(x) + tn x (x) otherwise 


(9.13) 


is well defined and is a C p 1 diffeomorphism on an open neighborhood of dX in M 3 (see Figure [Tt]) . 


Moreover, Ext maps dX x {0} x {0} to the boundary dX of A and for 0 < 77 < 770 if 

X~ v = X \ Ext(9A x] - 77 , 0] x {0}) 

X~ p is a compact C p_1 submanifold of A whereas 

X +p = A U Ext(<9Ax]0, + 77 ] x {0}) (9.14) 

is a C p_1 compact 2 dimensional submanifold of M 3 extending A along its boundary. 

Next theorem proves an approximation result by using C' p_ 1 (A)-functions. 
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Ext(x, s, i) 



(a) s < 0 (b) s > 0 


Figure 14: The map Ext 

Theorem 9.5. Let X be a C p (p ^ 2) compact oriented 2D submanifold of M 3 with non empty 
boundary denoted dX. Let / E BV(X, M) and let e > 0. Then there exists r\ > 0 and / E 
CP- 1 (X +p ,R d ) such that 


lx 


\f-f\dn 2 + \D x f\(X)-\D x f\(X) 


£ 


and 

[ I/I + \y x +vf\dn 2 ^ e. 

Jx+vXX-n 

Proof. We have first the existence of /' E C' p_1 (Xo,M) such that 

/ |/-/'|^ 2 + IDx/K^-IDx/'KX) 

Jx 

Moreover, there exists 0 < rf < t]q such that 

[ \f'\ + |Vx/ 7 | dT-L 2 < e 

JX\X-n 


^ e. 


Now, consider for r/ = r]'/3 the function f" : X +p —> M defined as 


/"(*) 


/(a) if z E 

/(Ext(x, —2r/ — s, 0)) if z E X +T? \ X~ v and where z = Ext(x, s, 0) 


Let <5 > 0, we check easily that for rf small enough 


Ix+^x-n 


I f 


+ |V*/"| dn 2 ^(l+5)x2 f \f'\ + \V x f 

Jxxx-^i 


dU 2 


since considering the change of variable induced by Ext we notice that (x,s) —>• ir x (x + sv(x)) 

has a determinant converging to 1 when (x, s ) convergence to a point on dX x {0} with s < 0. 

Moreover, we check easily that f"* _ E C P ~ 1 (X~ V ) and f"< _ E C p ~ 1 (X +r i \ X~p) where 

' x v 'x+v\x~v 

the intersection X~ p n X + p \ X~p = dX~ v is a C p ~ l one dimensional submanifold. Applying a 
smoothing in the vicinity of dX~ p we can obtain / E C 1 (A +?? ,M) such that / = f" on X~ v and 
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fx+v\x-v \ f\ + |Vx/| dH 2 ^ f x+ , 7 \ x - v | f"\ + \Vxf"\ dH 2 + £ so that (choosing 5 = 1 and rf small 
enough) 

|/| + |Vx/| cm 2 < 5e 


IX+v\X~v 


We get immediately 


/x 


|/-/|dW 2 + \D x f\(X)-\D x f\(X) < / |/-/'|dW 2 + px/IPO-II/v/'lpO 


IX 


and the result is proven. 


5e ^ 6e 


□ 


Remark 9.6 (L 2 and H 1 norms). Theorems 9.4 and 9.5 can also be stated by using the L 2 and 
H 1 norms instead of the 5b-one. The proofs in this case are simpler and follow the same argument 
than the previous ones. 
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